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1. Introduction and Summary 


In 2009 Petr Horava formulated new proposal of quantum theory of gravity (now known as 
Horava-Lifshitz gravity (HL gravity)) that is power counting renormalizable ill that 
is also expected that it reduces do General Relativity in the infrared (IR) limit 1 . The 
HL gravity is based on an idea that the Lorentz symmetry is restored in IR limit of given 
theory while it is absent in its high energy regime. For that reason Horava considered 
systems whose scaling at short distances exhibits a strong anisotropy between space and 
time, 

x ; = Zx , t' = l z t . (1-1) 

In order to have power counting renormalizable theory we have to demand that z > 3 in 
(3 + 1) dimensional space-time. It turns out however that the symmetry group of given 
theory is reduced from the full diffeomorphism invariance of General Relativity to the 
foliation preserving diffeomorphism 


x' 1 = x l + C(t, x) , t' = t + /(f) . 


( 1 . 2 ) 


Due to the fact that the diffeomorphism is restricted ( |1.2| ) one more degree of freedom 
appears that is a spin—0 graviton. The existence of this mode could have very significant 
consequences either for the consistency of given theory or for the phenomenological appli¬ 
cations of HL gravity. For that reason it would be desirable to formulate HL gravity where 
the number of the physical degrees of freedom is the same as in case of General Relativity. 
Such a proposal was formulated by Horava and Malby-Thompson in 115] in the context of 
the projectable HL gravity 2 . Their construction is based on an extension of the foliation 
preserving diffeomorphism in such a way that the theory is invariant under additional local 
U( 1) symmetry. The resulting theory is known as non-relativistic covariant theory of grav¬ 
ity. It was shown in [15, [l(| that the presence of this new symmetry implies that the spin-0 
graviton becomes non-propagating and the spectrum of the linear fluctuations around the 
background solution coincides with the fluctuation spectrum of General Relativity. 


1 For review and extensive list of references, see BI00- 

2 See also (||] and [Q 
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It is also well known that General Relativity contains large number of symmetries. 
Fixing all these symmetries we find that there are only two physical degrees of freedom 
left. Then we can ask the question whether it is possible to formulate the action for 
these physical degrees of freedom that is not based on the principle of covariance of the 
action under general diffeomorphism. The construction of such an action was proposed 
recently in two very interesting papers [^, 17]. The basic idea presented there was to 
perform the conformal traceless decomposition of the gravitational field [p] so that we 
have one degrees of freedom corresponding to the scale factor of the metric while we 
have five degrees of freedom of the metric that is restricted to have unit determinant. 
Then it was shown in Jl8|] that by gauge fixing of the Hamiltonian constraint one can 
eliminate the scale factor together with the conjugate momenta. As a result we obtain the 
action for five degrees of freedom that is invariant under spatial diffeomorphism where now 
Hamiltonian is determined by the solving of the Hamiltonian constraint of the General 
Relativity for the momentum conjugate to the scale factor. This analysis was further 
generalized in a very nice paper in [17| where the starting point was the action for the five 
physical degrees of freedom where it is required that given theory is invariant under spatial 
diffeomorphism. In other words we demand that the constraints that are generators of the 
spatial diffeomorphism are the first class constraints. We also have to require that these 
generators are preserved during the time evolution of the system. Then the requirement 
of the closure of the algebra of the Poisson brackets of these constraints together with the 
requirement of their time preservation determines the form of the Hamiltonian and the 
form of these constraints. When it is presumed that these constraints depend on partial 
derivatives of g l3 trough the scalar curvature we find that the original General Relativity 
action is reproduced. 

The goal of this paper is to formulate HL gravity for the gravitational physical degrees 
of freedom only in the similar way as in |l8|. To do this we start with another version of 
HL gravity that has the correct number of physical degrees of freedom and which is known 
as Lagrange multiplier modified HL gravity []7j. This model is based on the formulation of 
the HL gravity with reduced symmetry group known as restricted-foliation-preserving Diff 
(RFDiff) HL gravity 00- This is the theory that is invariant under following symmetries 


t' = t + 5t ,5t = const , x n = x 1 + C*(x, 


(1.3) 


The characteristic property of Lagrange multiplier modified HL gravity is an absence of the 


Hamiltonian constraint 113] and also presence of the additional constraint which changes 
the constraint structure of given theory so that the number of physical degrees of freedom 
is the same as in the case of General Relativity. Then in order to separate physical degrees 
of freedom of HL gravity we perform conformal traceless decomposition of the gravita¬ 
tional field, following || |8|. In this procedure we introduce new additional scalar field with 
additional symmetry so that the number of physical degrees of freedom is the same. Per¬ 
forming Hamiltonian analysis we also identify two second class constraints that, together 
with the gauge fixing scaling symmetry allow us to find Hamiltonian for the physical de¬ 
grees of freedom, at least in principle. These physical degrees of freedom are metric with 
unit determinant and conjugate traceless momenta so that the number of physical degrees 
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of freedom is the same as in the case of General Relativity. On the other hand there are 
also important differences. Since this theory arises from the theory with the complicated 
second class constraints we find that there is a very complicated symplectic structure on 
the phase space of the physical degrees of freedom. Secondly, even if we can claim that 
these second class constraints can be solved in principle we find that their solutions have 
the form of the non-local perturbative expansions. In other words it is hard to see how 
such a theory could be useful for some practical computations or even for its path integral 
formulation. 

However we mean that the analysis performed here suggests very interesting direc¬ 
tion in further research. The starting point would be the general form of the action for 
the physical degrees of freedom as was analyzed in [17] where we now presume that the 
additional term in the diffeomorphism constraint depends either on higher order of scalar 
curvature as for example R VJ R lJ or it depends on Rij non-locally. Then we should proceed 
as in [17] where we demand that the Poisson brackets of the spatial diffeomorphism con¬ 
straints close on the constraint surface. Then from the requirement of the preservation of 
these constraints during the time evolution of the system we could determine corresponding 
Hamiltonian density. We hope to return to this problem in future. 


2. Brief Review of Lagrange Multiplier Modified HL Gravity 


We begin this section with the brief review of the Lagrange multiplier modified RFDiff 
invariant HL gravity, for more detailed treatment see 0. RFDiff invariant Horava-Lifshitz 
gravity was introduced in [12], see also |L|. In 0 this action was extended by introducing 
Lagrange multiplier term that ensures that the spatial curvature is constant. Explic¬ 
itly,Lagrange multiplier modified RFDiff HL gravity has the form 


1 


s = — dtd^Vh{k i3 g ljkl k kl - v(h) + g[R\A ), 


( 2 . 1 ) 


where g[R ] = R — D, where D is constant, A is Lagrange multiplier that transforms as 
scalar 

A'(t',x.') = A(t,x) ( 2 . 2 ) 

under (L3). Further, K{j introduced in (p.l|) is modified extrinsic curvature 


IUj = -(dthij - ViNj - VjNi) 


(2.3) 


that differs from the standard extrinsic curvature by absence of the lapse N(t). Further 
the generalized De Witt metric g ijld is defined as 


gijki = }_( h ik h ji + h il h jkj _ Xh ij h ki ^ ( 2 . 4 ) 

where A is a real constant that in case of General Relativity is equal to one. Finally 
V(h) is a general function of hjj and its covariant derivative. The analysis performed in 
|7j showed that this theory possesses the same number of physical degrees of freedom as 
General Relativity. For that reason we mean that this action is a good candidate for the 
conformal traceless decomposition of the gravitational field and possible identification of 
the physical degrees of freedom of HL gravity. 
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3. Conformal Traceless Decomposition 


The conformal-traceless decomposition of the gravitational field was firstly performed in 
H in its initial value problem 3 . In order to implement conformal-traceless decomposition 
we follow [S] and define and Kjj as 

hij — (p (jij i Kij — (j) Aij + — (j) C/ijT . (’1 ■ 1) 

We see that this definition is redundant since the multiple of the fields g t j , q i>, , r give 
the same physical metric h k j and modified extrinsic curvature Kjj. In fact, we see that the 
decomposition (|3.1| ) is invariant under the conformal transformation 

= ^ 4 (x,t)^(x,t) , </>'(x,t) = n _ 1 (x,i)</>(x,i) , 

A-'ijfct) = n~ 2 {x,t)A i:j (x,t), r'(x,t) = r(x, t) . 

(3.2) 


We also see that (3J) is invariant under following transformation 


r '(x, t) = r(x,t) + C(x,t) , A l ij (yL,t) = A i:j (x,t) - ^C(x,t)^ 6 fty(x, t) . (3.3) 

Clearly the gauge fixing of these symmetries we can eliminate r and </>. 

In order to perform the Hamiltonian analysis of the conformal decomposition of the 
action m we firstly rewrite the action ®> to its Hamiltonian form. To do this we 
introduce the conjugate momenta 


pij _ 


<5,S 

Sdfhij 


= \VhG ijkl K kl , Pi = 


5S 

5d t N* 


= 0 , Pa = 


5S 

5d t A 


Then we easily determine corresponding Hamiltonian 


H = 


d'’’x(()ihi j P 1 ' - C) = J d 3 x(W T + WU'i) 


where 


u' T = ^P ij GijkiP kl + VgV(h) - VhAG(R ), n[ = -2h tJ v k p jk 


(3.4) 


(3.5) 


(3.6) 


Using the Hamiltonian and the corresponding canonical variables we write the action ( 0 ) 


as 


S = J dtL = J dtd 3 x{P ij d t hij - H) = J dtd 3 x(P ij d t hij - NU' T - iTHfi . (3.7) 

Then we insert the decomposition Q into the definition of the canonical momenta P l i 

P ij = 4 Va{^G ljkl A kl + \^rG ijkl g k i) , (3.8) 

O 


3 For review and extensive list of references, see |ld|]. 


- 4 - 







where the metric Q l i kl is defined as 


Q ijkl = g ik g> 1 + g il gi k ) - Xg ij g kl , Q^ kl = ^~ 8 g^ kl . 


Note that g i i kl has the i 


inverse 


Qijkl — 2 (9ik9jl d" 9il9jk ) j 9ij9kl ? Qijkl — 0 Qijkl • 

Using ( |3~8|) and fj we rewrite P i] df hij into the form 

P lJ d t h l3 = (jL^ggv kl A kl + - 3X)rgA d t9ji + 

+ ( ■^y/g^r l A k ig kl {l - 3A) + -^(1 - 3A)</> 5 4 d t <j> 

\ Hi kv J 


(3.9) 


(3.10) 


(3.11) 


We see that it is natural to identify the expression in the parenthesis with momentum 7r* J 
conjugate to gij and p<p conjugate to (j) respectively 

7r ij = ^g ijkl A ki + ^§(1 - llXyfrgV , 


3k 2 

V4> = 4\/5^ lA d.9 Jl (l - 3A) + 4r(l - 3X)cf> 5 T 

K> hi 


Then using ( |3 .1 2| ) we obtain following primary constraint 

Po ■ V<f>(t> ~ Xpigji = 0 . 


(3.12) 


(3.13) 


As we will see below this is the constraint that generates conformal transformation of the 
dynamical fields. Further, using (3.12) we find the relation between P lJ and 7r y in the form 


P ij = . 

Then we find that the kinetic term in "Ht takes the form 

- 2 K 2 !^ -6 


jj=p l] g ijk iP kl = '-^=-^g ijkl x 


kl 


(3.14) 


(3.15) 


As the next step we introduce the decomposition ( |3.1|) into the contribution f d 8 ~x.N i 'H' i . 
Using the relation between Levi-Civita connections evaluated with the metric components 
hij and g%j 

Y%{h) = T%{g) + 2i(^J + d^5 k - d^g kl 9lj ) (3.16) 

and also if we define n* through the relation N t = (j) 4 n t we obtain 


J d 3 xA PUi = J (P^ri U" , 


ini" 


(3.17) 
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where 


n'l = - 2 g ik DjTr jk + A^d^gwr kl , ( 3 . 18 ) 

where the covariant derivative Di is defined using the Levi-Civita connection T k -{g). Ob¬ 
serve that with the help of the constraint Ep we can write the constraint H'! as 

R" = -2g ik D^ k + d t <p P4> - = n, ~ A^dt^D ( 3 . 19 ) 

so that we see that it is natural to identify Ri as an independent constraint. In fact, we 
will see that the smeared form of this constraint generates the spatial diffeomorphism. 

Finally we should proceed to the analysis of the spatial curvature and generally the 
whole potential term V. Note that this is the function of the covariant derivative, R and 
Rij. Using the following formulas 

Rij[h] = RiM + ^D iC t>D j 4>-‘^D i D j ^-2 9 -^D k [g kl D l cl>\ - ^ gij D k fa kl Drf , 

R[h] = r 4 [R[g\ - , 

( 3 . 20 ) 

Then using also the relation between Levi-Civita connections evaluated on h and g we find 
that the potential term is generally function of cj) and g whose explicit form is not needed 
here. As a result we find the action in the form 

S = J dtd?x.(TT lj d t gij + Ptfidtcj) - n l Ri - - n lj <3 ijk iir kl - y/g^V^, h ) + 

+ y/g<f> 6 AS(<f>- A R\9] - 4 9 ij DiDj<i >) - AE D ) , 

( 3 . 21 ) 


where we included the primary constraint E o multiplied by the Lagrange multiplier A. 

Now we can proceed to the Hamiltonian analysis of the conformal decomposition of the 
gravitational field given by the action ( 3.21|) . Clearly we have following primary constraints 


TTi 


0 , « 0 , Y, d « 0 


(3.22) 


where n Tj and 7 t_ 4 are momenta conjugate to ri 1 and A with following non-zero Poisson 
brackets 

{n*(x),7Tj(y)} = 5j6{x - y) , (Al(x), vr^(y)} = <5(x - y) . (3.23) 

Further, the preservation of the primary constraints 7Tj and it a implies following secondary 
ones 

Ri ~ o , = ^y/g^Q « 0 . (3.24) 

Hi 

Now we should analyze the requirement of the preservation of the primary constraint Epp 
during the time evolution of the system. First of all the explicit calculations give 


{E D (x),5j i (y)} = 4 3 ij(x)tf(x-y) , 
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(3.25) 


{S Z) (x),7r y (y)} = —4vr y (x)(5(x - y) , 
{S D (x),</>(y)} = -</>(x)<5(x- y) , 
{S n (x),p^(y)} = <f>(x)5(x — y) 


using the canonical Poisson brackets 

{flr/( x )>^(y)} = + SlSj)S(x - y) , {0(x),jty(y)} = 5(x - y) . (3.26) 

It turns out that it is useful to introduce the smeared forms of the constraints 77 *, Xd 

Ts(AP) = J d' i xN , 'H l , D(M) = J d 3 xMS D , (3.27) 

where N l and M are smooth functions on R 3 . Then using (3.25) and also 

{So(x),T|(y)} = 25 k j d y i5(x - y) + 2 Sfd yj 6(x - y) - 2g kl (y)d y i5{x - y)^(y) 


we easily find that 
where 


{D(M),?4(y)}=0, 

Wrr = g ijk i7T kl + . 


Va 


(3.28) 

(3.29) 

(3.30) 


To proceed further we use following Poisson brackets 


{T s(N l ),g i: j(x)} = -N k d k g i:j (x) - diN k g kj (x) - g ik djN k {x) 
{T^(AT*),7r^(x)} = - d k (N k Tr ij )(x) + d k N i Tr kj (x) + Tr ik d k N\x ) 
{Ts(AP)>(x)} = -isrdittx), 

{T s (i\P),^(x)} = -d^N'p^ 


x 


and hence it is easy to see that 

{T 5 (I\P),£d(x)} = -irdiXoix) - d i N i 'E D (x) 


(3.31) 


(3.32) 


that together with ( 3.29| ) implies that ~ 0 is the first class constraint. 

Now we proceed to the analysis of the preservation of the secondary constraints 77* ~ 0 
and d>i ~ 0. Note that the total Hamiltonian takes the form 


H t = / d 3 x(74 + A£ d + n*77* + y PA + T/^i) , 


(3.33) 


where 7 is the Lagrange multiplier corresponding to the constraint p A while T/ is the 
Lagrange multiplier corresponding to the constraint $1 ~ 0. 
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In case of R- t we find following Poisson brackets 


r\ n 


(3.34) 


which implies that the smeared form of the diffeomorphism constraints takes the familiar 
form 

{T s (A^ i ), T s {M j )} = T s (N j d j M i - M j d 3 N l ) . (3.35) 

Further using (j3.31|) we easily find 


{T 5 (AP),7^(x)} = —d i N i 'H , T {x) — N i diH , T {ji) , 

{T silT),= —9jiV*$ 1 (x) - jrd&xW 


(3.36) 


that implies that Ri are the first class constraints that are preserved during the time 
evolution of the system. Finally we analyze the time evolution of the constraint ~ 0 . 
Using following formulas 

{i?(x),7r ?J '(y)} = -R lj (x)8(x - y) + D l D j 8(x - y) - g tj D k D k 8(x - y) , 

{ Ff,- (x), ir mn (y)} = \g kp [D i 6(6?% + %8?)6(x-y) + 

+ Dj(8™5™ + Sp8™)8(x - y) - D P (S?% + S?8?)8(x - y)] 

(3.37) 


we find that the time derivative of <f>i is equal to 
d t $ i = 

- - 0 x A , Rn) + 


+ 


^Va 

2 n 2 (j> 2 


3A - 1 

D k DM- 6 TT kl ] + X ) D fc U> fc [0-6 tt] - 

k nr 1 3A-1 ^ J 


K 2 A 

— 16—=- (it 12 --- g l i Tr)DiDj(j) + 

4> 5 y/g'' 3A — l y ’ 1 

16 4>k 2 „ , „ r ,_ R 8 k 2 4> 2A — 1 


+ 


Va 


Di(t>Dj[(j) b i r 


- 6 7rbl - 


Di4>g lJ Dj[4> 6 7r] = <P 2 , 


Va 3 A-i 


(3.38) 


where 4*2 is an additional constraint that has to be imposed on the system. Following 
119, 2C, 21] we include the constraint $2 into the definition of the total Hamiltonian that 
now has the form 


Hx 


J d 3 x(R' T - y/_gcf) 6 A((l)- 4 R[g\ - 80“ DiDrf - 0) + AS D + 


+ n l iii + 7 pa + r 7 4>i + r//<f> 2 ) 


(3.39) 
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Now we should again check the stability of all constraints. It is easy to see that the primary 
constraints together with Ts(N l ) are preserved while the time evolution of the constraint 
<hi « 0 is equal to 


<9^1 = j d 3 x( r // (x){$ 1 ,$ 2 (x)}) 

~ j d 3 xr / 7 (x){$ 1 ,$ 2 (x)} = o . 

As follows from the explicit form of the constraints d>i ]2 we have 

{<f>i(x),<h 2 (y)} / 0 . 


(3.40) 


(3.41) 


Then we find that the equation ( 3.40| ) gives T 11 = 0. In the same way the requirement of 
the preservation of the constraint <h 2 implies 


d t <&2 ~ / d x({4> 2 ,^t(x)} + T j (x){$ 2 ,4>i(x)}) = 0 . 


(3.42) 


Using the fact that {<h 2 ,%r(x)} / 0 and also the equation (3.41) we see that ( p.42| ) can 
be solved for T 1 . In fact, (|3.41 ) shows that dq and d> 2 are the second class constraints. We 
also see from the previous analysis that no additional constraints have to be imposed on 
the system. 

In order to find the action for the physical degrees of freedom we have to finally fix 
the gauge symmetry generated by Ti£>. To do this we introduce the gauge fixing function 


®G.F. = \[g- 1 • (3.43) 

It is easy to see that there is non-zero Poisson bracket between &g.f. and Ti£> so that they 
are the second class constraints. In summary we have following collection of the second 
class constraints 

4*1 = 0 , <3? 2 = 0 , = 0 , &g.f. = 0 • (3.44) 

The goal is to eliminate some degrees of freedom from these constraints, at least in principle. 
In fact, from dq, which is version of Lichnerowitz-York equation [^], we express <f> as 

^iv- 1 ^]-^) . (3.45) 

where V -1 is inverse operator to g l i D, t Dj. We can solve the equation above perturbatively 
around some constant 4> o- Further, from S/) we express p Finally, &g.f. reduces number 
of degrees of freedom in g to be equal to five and from <h 2 we express n as the function of 
remaining degrees of freedom. In summary, the physical degrees of freedom of Lagrange 
multiplier modified HL gravity are 

9ij , Vg = 1 , ,gijTr Jl = 0 . (3.46) 
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Note that there are three first class constraints Hi so that by gauge fixing these constraints 
we should eliminate remaining three degrees of freedom in gj 3 . In other words the physical 
content of given theory is the same as in General Relativity. 

From the previous analysis we see that 0 is a non-local function of R. The same 
situation also occurs in case of it so that in principle Hi has the form 


;Hi = -2^V^ + P 4V- i /(ff,vr))V^(V- i J R( 5 )) - ^(V 1 ,^) 


(3.47) 


By definition they are the first class constraints which however depends non-locally on 
the canonical variables. We would like to stress that this result can be considered as the 
generalization of the analysis performed in [17, 181 to the case of HL gravity. In more 


details, papers cited above were devoted to the construction of the action for the physical 
modes of the gravitational fields only that are the metric that obeys the condition ^fg = 1 
and the conjugate momenta n lJ that are traceless. In order to have the right number of 
physical degrees of freedom this action should be invariant under spatial diffeomorphism. 
Explicitly, the action studied there has the form 


S = J dtd 3 x(gijjt 1 ^ — ttr — N l Hi) , (3.48) 

where 

Hi = — k^ k ~ V iTTR , (3.49) 

where ttk is arbitrary function that has to be determined in such a way that Hi are the 
first class constraints. However when we restrict to the case when this function depends 
on the partial derivatives of gij through the scalar curvature R we find that the functions 
7 t/v and 7 tk are uniquely determined and leads to the General Relativity action. Note also 
that the symplectic structure used in given paper has the standard form 

{an,** 1 } = ^Sty + Sfl) - \ 9ij g kl . (3.50) 

In case of Lagrange multiplier modified HL gravity the situation is more involved. The 
symplectic structure is determined by the Poisson brackets of the second class constraints 
<J>i, <I >2 which is rather complicated. More precisely, let us denote all second class constraints 
as <3 ?,4 = (£_d, Tg.f. ■ 3>i, 4 > 2 ). Then the Poisson bracket between the constraints can be 
written as 

{4> j 4 (x),4> B (y)} = A AB (x,y) , (3.51) 

where the matrix A A b has following structure 


Aab(x, y) 


/o * o o\ 

* 0 0 * 

0 0 0 * 
\0 * * */ 


(3.52) 


where * denotes non-zero elements that depend on the phase space variables and their 
derivatives. Now it is easy to see that the Dirac brackets between the canonical variables 
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that are defined as 


{9ij(*),9kl(y)} D = ~ j dzdz' {gij(x), & A (z)} (A l ) AB (z,z') '),9kl{y)} , 

|vr^(x),7r“(y)|^ = - J dzdz {tt u (x), $a(z)} (A~ 1 ) ab (z, z') {$b(z') ; 7r w (y)| , 
{sij(x), 7r fcI (y)} D = {^(x),7r w (y)} - 
- J dzdz / {^(x),$ A (z)}(A _1 ) j4B (z,z , ){$s(z'),7r w (y)} 

(3.53) 


depend on the phase-space variables. Secondly, the resulting action is non-local due the 
presence of the inverse operator V” 1 . 

Even if the action for the physical degrees of freedom of Lagrange multiplier modified 
HL gravity is rather involved we mean that the results derived here should be considered 
as the starting point for further research of HL gravity when we generalize the analysis 
performed in [17, in several different ways. We could start with the action for the 
physical modes with the symplectic structure given by the equation ( 3.5C| ) and presume 
that 7 tk depends on R non-locally and try to determine the original action. Another 
possibility is to consider the case when ttk in ( |3.49 ) depends on R. l3 R l] and covariant 
derivatives of Rij. We hope to return to these problems in future. 
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